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(54) Monitoring the chemical composlion of a fluid 

(57) Apparatus for monitoring the chemical composition of a fluid eg liquid 14 having at least two components comprises 
sensors such as temperature probe 16. pressure dip tubes 20 for density measurement, and conductivity probe 18 for 
measurir^ at least two physical parameters dependent on the relative proportions of the components. The electrical 
measurements are processed by a micro processor using a series of polynomial equations relating the parameters with 
concentration and which are solved simultaneously by an iterative technique to derive values of the concentrations of the 
different components in the fluid. Control means may receive feedback signals so that any drift from a desired cone ntration 
can be corrected. The apparatus may be used in a plant for dissolving -uranium ore in nitric acid, the components being acid 
and Wjater. or In reprTOessI^^^ . 
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Monitoring the chemical composition of a fluid 

This invention relates to the monitoring of the 
chemical composition of a fluid, eg a liquor, in which at 
least two components are present. The liquor may for 
example comprise an acid solution of a uranium salt as 
obtained in the course of reprocessing irradiated nuclear 
fuels in which case the role of the invention may lie in 
the monitoring of the acid and uranium concentrations so 
that appropriate control actions can be taken to optimise 
the process. 

One method of determining acidity and concentration 
values is to use a titration technique but this suffers 
from the drawback that samples of the liquor have to be 
taken periodically and processing of the samples using 
titration techniques inevitably yields results some 
considerable time after taking the samples. 

British Patent No 950198 discloses an alternative 
approach in which concentration measurements are made by 
measuring certain "global" physical parameters of the 
liquor, namely conductivity and specific gravity, these 
being parameters which are governed by the relative 
concentrations of the acid and uranium components of the 
liquor. The apparatus described in the prior Patent 
relied on the fact that, while the relationships between 
conductivity, specific gravity, acid and uranium 
concentration are not linear, linear approximations can 
be employed when working within a limited range centred . 
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n predetermined concentration values and, in practice, a 
differently adjusted pair of resolver networks had to be 
employed for each working range* The known apparatus was 
also limited to a narrow range of temperatures (25 

5 5'C)r temperature compensation being affected by scaling 

the electrical analogue signals corresponding to 
conductivity and specific gravity respectively in 
electrical networks incorporating resistance 
thermometers. There is also the major problem of 

10. interaction between the components in the system, which 
is dealt with only in a very limited way by the resolver 
networks of Patent No 950198. Such interaction results 
in responses that are not additive; for example the 
conductivity of a solution containing a given uranium and 

15 free acid concentration, at a given temperature, is not 
the sum of the conductivities of the two solutions 
containing the same concentrations' of ~u~ralfi~um 'an~d~"f r e"e 
acid separately. 

The object of the present invention is to provide 

20 improved monitoring of the chemical composition of a 
fluid over extensive ranges of concentration and 
temperature where the relationships may be non-linear 
with complex interactions between the component 
concentrations and/or temperature. 

25 According to one aspect of the present invention 

ther is provided apparatus for monitoring the chemical 
composition of a fluid, comprising sensor means for 
providing respective electrical measures of at least two 
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physical parameters of the fluid, which parameters are 
dependent on the relative proportions of the components 
present in the systemr means for converting said 
electrical measures into digital form and digital 
computing means for utilising the resulting digital 
values to derive, by an iterative technique, values of 

g M2 • • • satisfying simultaneous equations of the 
form f^ ^ 9n (M-i, M2 • • •) -Pu where fjj = 0 

is the digital value of the Nth parameter and g^^ 
are second or higher order polynomial functions relating 
the concentrations , M2 • • • of the components 
in the fluid to the respective parameter value and 
where the functions g^ contain cross product terms, 
representing interactions, involving the concentrations 

, f^2 • • • 

According to a second aspect of the invention there 
is provided fluid processing plant in which the fluid 
comprises at least two components the concentrations of 
which are to be controlled, said plant including: me'ans 
for controlling the Introduction of at least one of the 
components so that the concentration of at least one of 
the components is thereby varied, sensor means for 
providing respective electrical measures of at least two 
physical parameters of the fluid, which parameters are 
dependent on the concentrations of the components present 
in the system, means for converting said electrical 
measures into digital form and 



digi-kal computing means for utilising the resulting 
digital values to derive^ by an iterative technique, 
values of H^, M2 > • • satisfying simultaneous 
equations of the form fj^ « gjj (M-i, ^2 • • "^N 

5 where « O, Pn is the digital value of the Kth 

parameter and g|| are second or higher order polynomial 
functions relating the concentrations M'l , M2 • • • 
of the components in the fluid to the respective 
parameter value Pjj and where the functions g^ contain 

10. cross product terms, representing interactions, involving 
the concentrations , M2 • • 

Preferably the plant or apparatus further includes 
means for sensing the temperature of the fluid and 
providing an electrical measure of the temperature (t) 

15 for conversion into digital form by the converting means, 

tlL© _f ^.i^ctiojis gjj in _this case also including t as a 

variable. 

In one example of the invention involving a two 
component liquor, eg containing nitric acid and uranyl 
20 nitrate, the equations can be of the form: 
f^ = Cq + C1M1 + C2M2 + C3t + + 
C5M2^ + C6t2 + C7H1M2 CgMit + 
C9M2t - Pi . . . ( 1 ) 

^2 ^o + "^l"! + ^^2^2 + "^3^ + ^4^1^ + 
25 K5M22 + J^e^^ + K:7MiM2 + KgH^t 

+ K9M2t - P2 . . . (2 ) 
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where the coefficients - Co dnd Kq — Kg are 

constants which may be derived using multiple polynomial 
regression by measuring the customary dependent variables 
and P2 experimentally from a series of liquor 
5 samples covering a wide range of known concentrations 
(Hi, M2 ) of the components and known temperatures, 

, M2 and t being the customary independent 
variables. Although and f2 are expressed above as 
second order polynomials, they may be expressed as third 
10 or higher order polynomials depending on the system, eg 
second order polynomials have proved adequate for a 
nitric acid/uranyl nitrate system whereas third order 
polynomials have been found to be more appropriate for 
certain concentration ranges of a KF/KOH system. 
15 Thus r having determined values for the co-efficients 

CQ-C9 and K^-Kg, substitution of the digital 
values P^i , P2 and t derived for a liquor with unknown 
component concentrations allows the two simultaneous 
equations to be solved iteratively for and M2- 

20 The iterative technique employed may involve the 

known Newton Raphson process but modified for the 
solution of simultaneous equations. 

For a two component system, the digital computing 
means may include: 

25 
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a* digital storage means for storing the values of the 
coefficients of the equations t^, ±2 which will be 
predetermined (by experiment) for a given system such as 
nitric acid/uranyl nitrate or KP/KOH; 
5 b. means for registering the outputs , P2 derived 
from the sensor means and, if applicable, the 
temperature prevailing in the system; 

c. means for holding provisional values Xi , X2 for 
Ml, M2; 

10 d. means for combining the provisional values Xi , 

X2 (and, if applicable, the measured temperature) 
with said coefficients in accordance with the 
respective polynomial equations f 1 , f 2 and 
calculating the resulting values of fi and f2r viz 

15 Fi and 

e;^ m~e"ans~f or-^combining "the -provisional -values xy," X2 

(and, if applicable, the measured temperature) with said 
coefficients in accordance with the respective first 
derivatives of fi, ±2 with respect to Mi, M2 in 

20 turn and calculating the resulting values Fi 1 , 

F12 and F2ir F22^ 
f • means for deriving from said values correction 
terms Ti , T2 of the form 

Ti « (F2*Fl2-Fi*F22)/(Fl 1*F22-F21*F12) 
25 T2 « (Fi*F2l-F2*Fn >/(Fll*F22-F21*Fl2)' 
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g* means for applying the correction terms and T2 to 

Xi and X2 respectively; and 

h. means for comparing the newly corrected Xi , X2 

values with the previously held X^ , ^2 values and 
depending on whether or not the result of comparison is 
favourable either updating the holding means with the 
corrected values and initiating a further cycle of 
iterative operation by means (d)-(h) or providing an 
output to indicate completion of the iterative process 
with the last corrected X^, X2 values representing 
the intera ti vely-deri ved final values for the 
concentrations , ^2* 

Thusr means (g) may be operable to form the sums 
Xi + and X2 + T2 whilst means (h) may effect 
comparison by determining the extent to which 
(Xi + Ti)/Xi and (X2 + ^2)7X2 deviate from 
unity; for example, if these ratios lie within the range 
0.9999 to 1#0001 the comparison may be designated 
favourable whereas ratios outside this range may be 
designated unfavourable • 

The values of Ml , M2 • • • may be employed in 
completing a feedback loop to allow the 

concentrations of the components to be maintained within 
predetermined limits* Thus » the controlling means is 
preferably supplied with feedback signals representing 
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the values M1 , M2 «. • • so that any drift away from the 
desired concentration levels in the system is corrected 
for by appropriate operation of the controlling means. 

In an alternative arrangement where feedback control 

5 is not required, means may be provided for activating 

some form of alarm or modifying operation of the plant 
(eg by closure or opening of valves or such like) in the 
event that the concentrations of the components 

deviate from preselected working ranges. 

10 The "global" parameters sensed by the sensor means 

may be density and conductivity but other parameters r 
such as refractive index, surface tension, viscosity, 
speed of sound, light absorbence/transmissivity output 
counts from a gamma absorptiome ter , XRF spectrometer etc, 

15 may be employed depending on the fluid. As a general 
ruTe ,~ for good results , o^f the" par Mi^"te"r"s~itreas~ulred~at 
least one should have a response to variations in 
component concentration "which is significantly different 
for the components. For example, in the case of a fluid 

20 containing nitric acid/uranyl nitrate, density variations 
are responsive predominantly to changes in uranyl 
nitrate concentration whereas nitric acid concentration 
changes tend to predominate in the case of conductivity 
variations • 

25 To promote further understanding of the invention, 

examples of the invention will now be described with 
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reference to the accompanying drawings in which: 

Figure 1 indicates diagramma ti cally a multi*probe 
sensor arrangement* for sensing density, conductivity and 
temperature of a KOH/KF solution contained within a 
vessel forming part of a HF scrubber using KOH as 
scrubbing solution; and 

Figure 2 is a schematic diagram showing a chemical 
plant employing a monitoring apparatus as part of a 
feedback loop for regulating concentration of a liquor. 

As shown in Figure 1 , the probe assembly comprises a 
protective vented tubular sheath 10 extending from the 
vessel top 12 and below the liquor level 14, The sheath 
10 houses a solid state temperature probe 16, an 
e lec tr ode less conductivity probe 18 and pneumicator tubes 
20 forming part of a density-measuring instrument. The 
latter operates by pumping a suitable gas such as 
nitrogen through the dip tubes 20 whose lower ends are 
.vertically staggered by a known distance and measures the 
difference in gas pressure, the density of the liquor 
being calculated from this difference. 

The analogue outputs from the probes 16, 18 and the 
density meter are fed via interface circuitry including 
analogue/digital converter circuits to a microprocessor 
programmed to execute the mathematical routines that will 
be described below. These routines are based on the 
polynomial equations (1) and (2) set out hereinbefore. 
As mentioned previously^ the coefficients/constants 
involved in these equations must be derived 



1 0 

experimentally for a given system, for example by making 
up sample solutions containing a wide range of known 
<5oncentrations of the components (ie covering at least 
the range that can be expected in practice) and measuring 

5 the "global" parameters , ie density and conductivity in 
this case, for each sample over a range of temperatures 
(again covering at least the range that can be expected 
in practice). The experimental results obtained for 
Pi , P2 and the known concentrations M-| , M2 can 

10 then be processed using multiple regression analysis to 
obtain "best-fit" values for the coefficients C^, Ci 
. • • and Kq, K-| • . These values then apply for 
the given system and can be stored as an array in a 
digital memory associated with a microprocessor. 

15 In theory the C and K values can be derived by 

^ea:suri"ng~n~ so lu tion s -if- ther e- ar e— n— va lues -of — G— and— K 

in each of the 2 equations. Each of the n solutions 
would have a "different KOH and/or KF concentration and a 
different temperature. 

20 In practice, more than n solutions are always 

measured and each solution-is measured at different 
temperatures. This is convenient practically and gives 
more precise determinations of C and K. <The smaller n 
the less precise the determinations of C and K and hence 

25 the less precise the subsequently measured component 
concentrations). 

Typically 20-40 solutions are measured, each at 5-20 
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different temperatures covering the temperature range the 
plant solution will encounter. This results in hundreds 
of equations, which are solved for the C and K values 
using multiple polynomial regression. Selecting the 
5 concentrations of the 20-40 solutions is done as follows. 
{ The concentration range of each component occurring in 
the plant liquor is recorded and then extended somewhat 
at the lower and upper end. Each continuous range is 
then divided^ typically^ into 10-20 discrete 
10 concentrations. Each discrete concentration of one 

component is then linked either a) at random or b) in a 
systematic manner, depending on the application, to a 
discrete concentration of another component. Constraints 
are frequently put on this linking. These constraints 
15 are either a) a combination of component concentrations 
which is physically impossible or b) a combination which 
could not occur on the plant. 

For the HF scrubber, the 2 components are KOH and 
KF. The KOH absorbs HF which is converted to KP. 
20 Thus the KP concentration rises as the KOH 

concentration falls. For example, 28 solutions initially 
measured to determine the values of C and K were chosen 
systematically as follows, each solution being measured 
over the range 5-40*C at approximately 5*C increments. 
25 The conductivity and density of each solution was 

measured at the 8 different temperatures giving 224 
equations relating conductivity to KOH and KP 
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concentration and temperaturet and similarly for density. 
For a third order polynomial there are 19 values of C to 
determine (conductivity) and 19 values of K (density) • 
The solutions are always measured in a random order, to 
avoid systematic errors. 



SOLUTION 


qKOH/I 00ml 


<yKF/ 1 00ml 
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26 
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27 
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28 
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0 



C and K values are determined by multiple polynomial 
regression, using the method of least squares. In theory 
only 19 equations are required (for a third order 
polynomial system) to determine the 19 C values, and 
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similarly for the 1 9K values* This gives an 
over determined systemr ie more equations than 'necessary* 
and has the advantage o£ providing the uncertainty in the 
values of C and K, ie expressed as standard deviations* 
The method of least squares is typically used for solving 
overdeter mined systems • 

The microprocessor is connected to receive the 
outputs of the density, conductivity and temperature 
sensors via a suitable interface including 
ana logue -to-di gi tal conversion circuitry so that, in 
practice, the concentration levels in the liquor can be 
monitored by way of the measurements available for the 
density, conductivity and temperature. 

In practice for a KOH/KF system, third order 
polynomial equations have been used and the 
micropr,ocessor is programmed to execute the stepwise 
procedure outlined below. 

The following terms are used in- the equations that 
follow: 



y Conductivity (mS»cm''^), as measured 

w Density (g.cm"^), as measured 

t Temperature (deg C), as measured 

X KOH concentration (g KOH/IOOml) 

z KF concentration (g KF/IOOml) 

A Calculated conductivity 

B Calculated density 
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D Difference of calculated and measured 
conductivity 

B Difference of calculated and measured density 

P 1st differential of conductivity w.r.t. x 

5 R 1st differential of conductivity w.r.t. z 

Q 1st differential of density w.r.t. x 

S 1st differential of density w.r.t. z 

P Correction term for conductivity 

G Correction term for density 

10 El Error value for conductivity 

E2 Error value for density 



The above information is stored in memory arrays 
associated with the microprocessor to allow the 
intermediates to be printed at the end of a calculation. 

15 The coefficients for the conductivity and density 

equations are h^ld "in arralrs~C (D and" KM )- respectively-- - 
Initially, arbitrary values are selected for x and z. 
The temperature t is conveniently the difference between 
a centred value for the temperature (ie the median of the 

20 operating temperature* range of the plant) and the liquor 
temperature as measured. „ 

Step 1 - The values for x, z (initially arbitrary) and t 
are fed into the following equations to arrive 
at a calculated value of the conductivity (A) 
25 and density (B). 
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A « CiX+ C22 + C3t 

+ C4x2 + C5r2 f Cfit2 + C7XS5 + C'nxt + Cgsst 

+ C^5xz2 + C^gxt2 + C^yz2t + CjgZt^ + C,gX2t 

B « K-ix + + K3t 

+ K4x2 + K5Z2 + Kgt^ + Kyxz + Ksxt + KgZt 

+ K|qx3 + K^|z3 + K^2*^*^ + ^^13X^2 + Kj^x^t 

+ K^5Xz2 + K|5Xt2 + K^yZ^t + K|3Zt2 + K^gXZt 

Step 2 - The difference between the derived and 
experimental conductivity and density 
measurements is calculated. 

D = A - y [ =r fy (XrZ) 1 Conductivity 
E « B - w [ = (x,z) ) Density 

The y and w equations contain the unknowns x and z, and 
are solved for x and z using a Kewton Raphson iterative 
process adapted for th^s solution of simultaneous 
equations • 

Step 3 - The equation for A is differentiated first 
w.r.t. 'x' and then w.r.t. *z* to give the 
respective first differentials P and R. 
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P = C| 

+ 2C4X + C7Z + Cgt 

+ -3C,0*^ + 2C13XZ + 2C^4Xt 4- C^jzS + C^gt2 + c^gZt 
5 R = C2 

+ 2C5Z + C7X + Cgt 

+ 3C|^22 + C,3x2 + 2C15XZ + 2C,yZt + C^gt^ + C^gXt 



Step 4 - The equation for B is differentiated first 
10 w.r.t* 'x' and then w.r.t. 'z^ to give the 

respective first differentials Q and S. 

Q = Ki 

+ 2K4X + K7Z + Ket 

4 5- - -+- -3K,.o_x2_ i^2K, 3 XZ j^^ + K^^t^ + K^^xt 

S = K2 

+ 2K5Z + K7X + Kgt 
20 + 3K, ,z2 + K,3X2 + 2K^ 5 xz + 2K^ ^ z t + g 1 2 + ^ xt 

Step 5 - The correction term for the KOH and KF 
concentrations are calculated. 

P = (E*R - D*S) / (P*S - Q*R) KOH concentration 
25 G = (D*Q - E*P) / {P*S - Q*R) KF concentration 

Step 6 - The correction terms are now added to the 
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previous arbitrary values of the KOH and KF 
concentrations • 

= Xi + F 

« Zi + G 

Step 6 - The error between consecutive calculations is 
calcula ted« 
El = Xi+i /Xi 
= Zi+i /zi 



Steps 1 to 6 are then repeated each time using the 
corrected values for x and z until El and E2 lie between 
for example 0.9999 and I.OOOI. The values for x and z 
then prevailing are then taken as the KOH and KF 

15 concentrations existing in the liquor at that time and at 
the particular temperature measured. 

In this manner, the KOH and KF concentrations can be 
continuously monitored and control of the system may be 
modified according to requirements. For example, the rate 

20 of supply of KOH to the scrubber may be controlled in 

dependence upon the concentration levels calculated so as 
to maintain the KOH/KF concentration levels 
within predetermined limits. 

Referring now to Figure 2, this shows 

25 diagrammatically plant for dissolving uranium ore in 

nitric acid, the ore being introduced into vessel 30 at a 
constant feed rate by means of a screw feeder 32 and the 
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acid being introduced via supply line 34 along with water 
fed via supply line 36. The contents of the vessel are 
agitated by means of stirrer 38. The supply lines 34, 36 
include motorised flow-setting valves 40, 42 for 

5 regulation of supply of acid and water to the vessel 30 
in order to maintain the concentrations of uranium and 
nitric acid within preset limits. 

To control the valves 40, 42, the concentration of 
the liquor is monitored by means of a 

10 microprocessor-based analysis unit 44 which receives 

input signals from a probe assembly 46 which may be of 
the form shown in Figure 1. Thus, the probe assembly 46 
provides electrical signals representing the density, 
conductivity and temperature prevailing in the liquor. 

15 The signals are processed by the analysis unit 44 in the 

manner" des"cr iTa^d^above— to -pro^v^^^ -two analogue output 

signals which represent respectively the uranium and acid 
concentrations of the liquor and are applied via inputs 
48, 50 to a programmable controller 52 (such as a Model 

20 6433 Programmable Signal Processor marketed by Turnbull 
Control Systems Limited). The controller 52 is 
programmed to process the input signals to derive 3-term 
output control signals for operating the valves 40, 42 so 
as to maintain the flow rates of acid and water at values 

25 such that the uranium and acid concentrations remain 
within limits programmed into the controller 52. 
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In addition to providing signals representing 
concentration levels to the controller 52, the analysis 
unit 44 is also connected to an alarm 56 and is 
programmed to actuate the alarm in certain cir cums tances r 
eg if a fault develops in the unit 44 or the probe 
assembly or if the concentration of either of the species 
being monitored falls outside preset limits. 

Although the invention is described above in 
relation to specific applications, ie KOH/KF 
concentration levels in a HF scrubber or uranium/acid 
concentration levels in a dissolver# it will be 
understood that the same hardware can be used on 
different plant if the microprocessor is reprogrammed 
appropriately • 

The invention is also applicable to systems in which 
there are three or more interacting components* Also the 
ambit of the invention is not restricted to liquors; gas, 
vapour and liquid systems are all considered to be within 
the scope of the invention. The parameters measured may 
be entirely different' physical entities such as density 
and conductivity but they may be different 
characteristics of the same physical entity, eg 
ultraviolet absorption at different wavelengths. 
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Claims 

1. Apparatus for monitoring the chemical composition of 
a fluid comprising sensor means for providing respective 
electrical measures of at least two physical parameters 

5 of the fluid r which parameters are dependent on the 

relative proportions of the components present in the 
fluid, means for converting said electrical measures into 
digital form, and digital computing means for utilising 
the resulting digital values to derive, by an iterative 

10 technique, values of M-| , M2 • • • satisfying 

simultaneous equations of the form f = gjj (W^ , 

M2 • • •) -Ph where f^ = 0, is the digital 

value of the Nth parameter and are second or higher 

order polynomial functions relating the concentrations of 

"~~T5 the -components. M:|., JM2_^ . • in the fluid to the 

respective parameter value Pjj and where the functions 
gjj contain cross product terms, representing 
interactions, involving the concentrations Mi , M2 • • 

20 2. Apparatus as claimed in Claim 1 further including 
means for sensing the temperature of the fluid and 
providing an electrical measure of the temperature (t) 
for conversion into digital form by the converting means, 
the equations gjj in this case also including t as a 

25 variable. 

3. Apparatus as claimed in Claim 1 or Claim 2 in which 
said digital computing means comprises; 
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digital storage means for storing the values of the 
coefficients of the polynomials fi, £2 v^^^^ 
will be predetermined (by experiment) for a given 
fluid such as nitric acid/uranyl nitrate or KF/KOH; 
means for registering the outputs P-| , P2 derived 
from the sensor means and, if applicabl-e, the 
temperature prevailing in the system; 
means for holding provisional values X-i r ^2 
Ml, M2; 

means for combining the provisional values X-i , 
X2 (andr if applicable, the measured temperature) 
with said coefficients in accordance with the 
respective polynomial functions f^, f2 and 
calculating the resulting values of f ^ , and f2f 
viz F-i and F2; 

means for combining the provisional values X-i , 
X2 (and, if applicable, the measured temperature) 
with said coefficients in accordance with the 
respective first derivatives of f 1 r f2 with 
respect to M-| , M2 in turn and calculating the 
resulting values F^ ^ , ^12' ^21 
^22? 

means for deriving from said values correction terms 

, T2 of the form 
Tt « (F2*F -Fi*F )/(F *F -F *F ) 
T2 = (F^*F -F2*F )/(F *r -F *F )y 
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g. ^^means for applying the correction factors and T2 

to Xi and X2 respectively; and 

h. means for comparing the newly corrected , X2 
values with the previously held X^ , X2 values and 

5 depending on whether or not the result of comparison is 

favourable either updating the holding means with the 
corrected values and initiating a further cycle of 
iterative operation by means (d)-(h) or providing an 
output to indicate completion of the iterative process 

10 with the last corrected Xj , X2 values representing 

the interatively-derived final values for the 
concentrations , M2 • 

4. In combination, apparatus as claimed in any one of 
Claims 1 to 3 and means operable to control the quantity of 

15 at least one of the components present in the fluid, the 

CO iTtrol" means being supplied- -with^f.eedbac^k__s.ignals from the 
monitoring apparatus so that any drift away from the desired 
concentration levels in the fluid can be detected by the 
monitoring apparatus and corrected for by appropriate 

20 operation of the control means. 

5. Apparatus for monitoring the chemical composition of a 
fluid substantially . as hereinbefore described with 
reference to, and as shown in, the accompanying drawing. 



25 



BNSOOCID: <GB 2217B48A_I_> 



23 

6. Fluid processing plant including apparatus a claimed 
in any preceding claims and including means for 
controlling the introduction of at least one of . the 
components so that the concentration of at least one of 
the components is therehy varied. 

3 0XD8 
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